/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      phaseProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

phases (RDX tnt gas);

RDX
{
    type detonating;
    reactants
    {
        thermoType
        {
            transport   const;
            thermo      eConst;
            equationOfState Murnaghan;
        }
        equationOfState
        {
            rho0            1160;
            n               7.4;
            kappa           3.9e11;
            Gamma           0.35;
            pRef            101298;
        }
        specie
        {
            molWeight       222.12;
        }
        transport
        {
            mu              0;              // Viscosity
            Pr              1;              // Prandtl number
        }
        thermodynamics
        {
            Cv              1130.78;           // Heat capacity
            Hf              0.0;
        }
    }
    products
    {
        thermoType
        {
            transport   const;
            thermo      eConst;
            equationOfState JWL;
        }
        equationOfState
        {
            A               2.9867e11;
            B               4.11706e9;
            R1              4.95;
            R2              1.15;
            C               7.206147e8;
            omega           0.35;
            rho0            1160;
        }
        specie
        {
            molWeight       55.0;
        }
        transport
        {
            mu              0;              // Viscosity
            Pr              1;              // Prandtl number
        }
        thermodynamics
        {
            Cv              1000;           // Heat capacity
            Hf              0.0;
        }
    }

    activationModel linear;
    initiation
    {
        E0      4.0e9;
        points  ((0 0.55 0));
        vDet    7850;
    }

    residualRho     1e-6;           // Minimum density of the phase
    residualAlpha   1e-6;          // Minimum volume fraction used for division
}

tnt
{
    type detonating;
    reactants
    {
        thermoType
        {
            transport   const;
            thermo      eConst;
            equationOfState BirchMurnaghan3;
        }
        equationOfState
        {
            rho0            1601;
            K0              9.6e9;
            K0Prime         6.6;
            Gamma           0.35;
            pRef            101298;
        }
        specie
        {
            molWeight       227.13;
        }
        transport
        {
            mu              0;              // Viscosity
            Pr              1;              // Prandtl number
        }
        thermodynamics
        {
            Cv              1095;           // Heat capacity
            Hf              0.0;
        }
    }
    products
    {
        thermoType
        {
            transport   const;
            thermo      eConst;
            equationOfState JWL;
        }
        equationOfState
        {
            rho0        1601;
            A           371.21e9;
            B           3.23e9;
            R1          4.15;
            R2          0.95;
            omega       0.3;
        }
        specie
        {
            molWeight       55.0;
        }
        transport
        {
            mu              0;              // Viscosity
            Pr              1;              // Prandtl number
        }
        thermodynamics
        {
            Cv              1000;           // Heat capacity
            Hf              0.0;
        }
    }

    activationModel pressureBased;
    initiation
    {
        E0              7.0e9;
        I               0;
        G1              3.5083e-7;
        c               1.0;
        d               0.0;
        y               1.3;
        minLambda1      0.0;

        G2              0.0;

        pMin            5e5; // Cutoff pressure for activation
    }

    residualRho     1e-6;           // Minimum density of the phase
    residualAlpha   1e-6;          // Minimum volume fraction used for division
}

gas
{
    type basic;
    thermoType
    {
        transport   const;
        thermo      eConst;
        equationOfState idealGas;
    }
    equationOfState
    {
        gamma           1.4;            // Heat capactiy ratio
    }
    specie
    {
        molWeight       28.97;
    }
    transport
    {
        mu              0;
        Pr              1;
    }
    thermodynamics
    {
        type        eConst;
        Cv          718;
        Hf          0;
    }

    residualRho     1e-6;
    residualAlpha   1e-6;
}


// ************************************************************************* //
